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Abstract 



The low energy dynamics of QCD is investigated with special attention paid 
to the matching between QCD and chiral perturbation theory(ChPT), and 
also to some useful algebraic chiral operator relations which survive even 
when we include chiral loop corrections. It then allows us to evaluate the 
hadronic matrix elements below the energy scale A x ~ 1 GeV. Based on the 
new analyzes, we present a consistent prediction for both direct CP-violating 
parameter e'/e and AI = 1/2 rule in the kaon decays. In the leading 1/N C 
approximation, the isospin amplitudes Aq and A2 are found to agree well 
with the data, and the direct CP-violating parameter e'/e is predicted to 
be large, which also confirms our early conclusion. Its numerical value is 
e'/e = 23.6l7 2 g 4 x 10~ 4 (ImXt/l-2 x 10~ 4 ) which is no longer sensitive to the 
strange quark mass due to the matching conditions. Taking into account a 
simultaneous consistent analysis on the isospin amplitudes Aq and A2, the 
ratio e'/e is in favor of the values e'/e = (20 ± 9) x 10 -4 . 

PACS numbers: 11.30.Er, 13.25.Es, 12.39.Fe 



1 



Typeset using REVTeX 



I. INTRODUCTION 



To make a consistent prediction for the direct CP- violating parameter e' / e caused by the 
Kobayashi-Maskawa CP-violating phase which can arise from the explicit CP violation 
in the standard model(SM) or originate from the spontaneous CP violation || in the simple 
extension of SM with two Higgs doublets (S2HDM) ||, it is thought to be necessary to 
understand simultaneously the longtime puzzle of the AI = 1/2 rule in kaon decays as they 
involve the long-distance evolution of common hadronic matrix elements. It is believed that 
the low energy dynamics of QCD shall play a crucial role for a consistent analysis. During 
the past few years, both theoretical and experimental efforts on direct CP violation in the 
kaon decays have been made important progresses. As a consequence, it has reached an 
agreement between the experimental results [§,^0] and the theoretical predictions within 
the framework of ChPT P,|T0f and chiral quark model [[II] . On the experimental side, 



two improved new experiments |||7j] with higher precision have reported results which are 
consistent each other and also agree with the early result [|J. On the theoretical side, 
there have been some developments which are mainly based on QCD of quarks and cut-off 
ChPT at low energies for mesons. The renormalization coefficients of all the relevant four 
quark operators, which characterize the short-distance effects of the effective Hamiltonian 
generated from renormalization of the weak interactions, has been computed and extended 



from the leading order [0,|TJ] to the next-to-leading order flqjlql QCD corrections. The 
results agree with each other. The long-distance effects have been evaluated from the ChPT 
inspired from 1/N C expansion [|17])0 U P ^° ^ ne chiral one- loop level [19|.|],|I(J as well as from 
chiral quark model Recently they have been recalculated within the same framework of 
ChPT but with a different calculating scheme [20 ^3|. The important issue concerned in all 
the calculations is the matching problem to QCD. In ChPT, it requires the matching between 
the short-distance operator evolution from QCD with infrared cut-off and long-distance 
operator evolution from ChPT with ultraviolet cut-off. Practically, the renormalization 
scale fi dependence of the Wilson coefficient functions Cj(/z) from perturbative QCD should 
cancel the one of the corresponding operators Qi{n) from non-perturbative contributions. 
In the chiral quark model the operator evolution leads to the results which are expected to 
valid only at a special value of the energy scale /i. Alternatively, in the ChPT approach, its 
attractive advantage is that chiral loops with an ultraviolet cut-off, denoted by M, introduce 
a scale-dependence for long-distance operator evolution. As a simple consideration, the 
ultraviolet cut-off M might naively be identified to the infrared cut-off fi to improve the 
matching. Consequently, both the AI = 1/2 rule and direct CP- violating parameter e'/e 
can be enhanced to be more consistent with the present experimental data. Nevertheless, in 
the existing treatments of the approach, there remain some open questions which need to be 
further clarified. Firstly, the momentum cut-off M in the long-distance operator evolution 
from meson loops can in general only be extended to the energy scale which must be smaller 
than the dynamical chiral symmetry breaking scale A/, i.e. M < Af ~ lGeV, whereas the 
short-distance operator evolution from perturbative QCD (by using renormalization group 
equation) requires that the energy scale should be above the confining scale, i.e., [i > 1 
GeV. Thus naively identifying the ultraviolet cut-off M in ChPT to the infrared cut-off \i in 
perturbative QCD may become inappropriate. Secondly, there appear some discrepancies 
between the ref. Jl9| and ref. |p0| , p2| for the matrix elements Q\ and Q2 even if the same 
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chiral Lagrangian has been used and the same loop diagrams have been considered. It is 
noticed that the discrepancies only occur in the coefficients of the quadratic terms of the cut- 
off energy scale M and in the constant terms. Such discrepancies mainly arise from different 



calculating schemes. In refs. [pLQ|,P|,[10|] , all the chiral one- loop contributions were considered 
to be summed up with the cut-off regularization, and the coupling constants are replaced by 
the renormalized ones. Such a treatment is the standard one as adopted in the quantum field 



theory. In the recent calculations [ETH23I] , the chiral one-loop diagrams have been separated 



into two classes, i.e., so-called factorized and nonfactorized diagrams. For the nonfactorized 
diagrams, a virtual momentum flow has been artificially added to the propagators. As a 
consequence, such an alternative treatment is equivalent to the change of cut-off energy scale. 
When taking the cut-off energy scale to be infinity as the case for a renormalized field theory 
like QCD and QED, the treatment has no effects. However, for a finite cut-off integral, the 
change of the variables of the integrand will result in different results. This is the main 
reason why two calculations led to different results in the quadratic terms of cut-off scale 
and in the constant terms. In fact, one can simply rescale the cut-off scale A c in refs. 



via 



K 2 C = 2M 2 /3 to obtain the results in ref. |19| for Qi and Q2- With a similar reason, for Q$ 
there also exist discrepancies between ref. [||[10| and ref. |23| for the quadratic and constant 
terms, however, two results cannot be simply related by the same rescaling factor as the one 
for Qi and Q2- The reason is that the results in ref. [f2^] were obtained only by evaluating 
part of the so-called nonfactorized diagrams. Notice that it may not be so clear to separate 
the factorized and nonfactorized diagrams for the loop corrections of the left-right structure 
four quark operators generated from the penguin diagrams in which the intermediate quarks 
form a closed loop with the gauge bosons. In fact, for the density x density operators, the 
so-called factorized diagrams do provide contributions to the anomalous dimension of the 
operators in QCD evolution. The //-dependence of such factorized contributions is exactly 
cancelled by the one of explicit quark mass factor appearing in the corresponding chiral 
operators, but not by the one of the corresponding factorized chiral loop, since the quark 
mass factor does not arise from the chiral loop contributions. Therefore, for the density 
x density operators, or more general for the operators with left-right structure, it is not 
necessary to have one to one correspondings between QCD loop and chiral loop due to the 
//-dependence of the low energy coupling constants in the chiral Lagrangian. Nevertheless, 
physics observables should be independent of the calculating schemes. 

Thus, clarifying the above two open questions comes to one of the main purposes of this 
paper. Our paper is organized as follows: in section 2, we describe the basic motivations 
for evaluating the long-distance contributions of the hadronic matrix elements within the 
framework of ChPT. Especially, a functional cut-off momentum is introduced for the purpose 
of matching, namely the cut-off momentum M is in general considered as the function of 
the QCD running scale /i, i.e., M = M(/i); in section 3 we explicitly write down the chiral 
representation of four qurak operators and emphasize some useful algebraic chiral operator 
relations; in section 4 we investigate the matching between QCD and ChPT, where the chiral 
operators are explicitly evaluated in the functional cut-off momentum scheme, two useful 
matching conditions will be obtained. Of interest, the strange quark mass is found to be 
fixed from the matching condition and algebraic chiral operator relation. The long-distance 
chiral operator evolution is carried out in section 5; in section 6 we present our numerical 
predictions for the direct CP-violating parameter e'/e and A J = 1/2 rule, they are found to 
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be remarkably consistent with the data. Our conclusions and remarks are made in the last 
section. 

II. BASIC MOTIVATION 

Our considerations are mainly based on the following basic points: 

• In the large N c limit but with the combination a s N c = «o being held fixed. The QCD 
loop corrections which are proportional to a s are then corresponding to a large N c 
expansion, a s ~ 1/N C |l7j. 

• Chiral symmetry is supposed to be broken dynamically due to attractive gauge interac- 
tions, namely the chiral condensates < qq > exist and lead to the Goldstone-like pseu- 
doscalar mesons 7r, K, rj. The chiral symmetry breaking scale Af is characterized by the 
condensate, Af ~ 47ry / — 2 < qq > jr ~ 1 GeV with r = m 2 o /rh (m = {m u + /2). 

• The chiral Lagrangian is considered to describe the low energy dynamics of QCD in 
large N c limit and is going to be treated as a cut-off effective field theory. The cut-off 
momentum M is expected to be below the chiral symmetry breaking scale Af. 

• The chiral meson loop contributions are characterized by the powers of p 2 /A 2 with 
Af = Airf . Here f 2 ~ —2 < qq > jr ~ N c is at the leading N c order and fixed by the 
7r decay coupling constant / ~ F^. Thus the chiral meson loop contributions are also 
corresponding to a large N c expansion of QCD, p 2 /A 2 ~ 1/N C ~ a s . Therefore both 
chiral loop and QCD loop contributions must be matched each other, at least in the 
sense of large N c limit. Thus the final physical results should be independent of the 
cut-off schemes. 

• The cut-off momentum M of loop integrals should not be naively identified to the 
renormalization scale fi appearing in the perturbative QCD in large N c limit. It is 
in general taken to be a function of /i, i.e., M = M(/x), which may be regarded as 
a functional cut-off momentum, its form is determined by the matching between the 
Wilson coefficients of QCD and hadronic matrix elements evaluated via ChPT. It is 
seen that the matching relates the chiral cut-off momentum to the strong coupling 
constant so that the results become scheme independent. 

From these points of view, the ChPT with functional cut-off momentum is going to 
be treated, in certain sense, as a low energy effective field theory of QCD in the large A^ c 
limit. With such a treatment, it is in general not necessary to distinguish the so-called 
factorized and non-factorized contributions since the renormalization of field theory should 
well cover both of their contributions automatically. In this paper, we will give up the 
calculating scheme of separating the factorized and nonfactorized contributions, and adopt 
the calculating scheme first proposed by Bardeen, Buras and Gerard [|l9) in the ChPT 
inspired by the 1/N C expansion, but with a functional cut-off momentum M(/x) instead of 
naively identifying the cut-off momentum to the QCD running scale [i. 
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III. CHIRAL REPRESENTATION AND ALGEBRAIC RELATIONS 



In the standard model, the AS = 1 low energy (// < m c ) effective Hamiltonian for 
calculating K — > im decay amplitudes can be written as 

Hffi x = %KY,c i ( J *)Q i {ji), (//<m c ) (1) 

j= i 

with Qj the four quark operators 



Qi 


= 4 s L Yd L ul^^ul , 


Q 2 


= 4 s l Yul u L lnd L , 


Q 3 




Qa 


= 4 J2 SlI^Ql QL7^d L , 


Q 5 


<? 


Qe 


= -8 SlQr QRd L , 
q 


Qi 


3 

q 1 


Qs 


3 

= -8 J2 n e q ^lQr qRd L 
q 1 



with q R L — ~(1 ± 7s)g. Where the sum goes over the light flavors (q = u, d, s) and e q is the 
charge of the corresponding light quarks. Q3, ■ ■ ■ ,Q& arise from strong penguin diagrams. 
They transform as (8^,1^) under SU(3)l x SU(3)r and solely contribute to AI = 1/2 
transitions. Note that only seven operators are independent as the linear relation Q4 = 
Q2 — Qi + Qz- Q7 an d Qs originate from electroweak penguin diagrams. are Wilson 

coefficient functions 



a(ji) = Zi(ji) + ryiiji) . (3) 

where r = —\ t /\ u with \ q = V*V qc i. The Wilson coefficient functions Zi(fi) and 
have been evaluated up to the next-to-leading order QCD corrections. The K — > Tin decay 
amplitudes Aj with isospin I are given by 

A^ 5 ' = {<K<K\H% s f\K) = % X^c^iQ^j (4) 

where 5i are the final state strong interaction phases. It is a hard task for calculating the 
hadronic matrix elements (Qi(n))i for fi < A x = 1 GeV which is at the order of chiral 
symmetry breaking scale. This is because perturbative QCD becomes unreliable in such a 
low energy scale. In this paper we adopt the ChPT with functional cut-off momentum to 
evaluate {Qi(fi))i when \x < A x . To do that, the procedure is as follows: one first represents 
the current x current or density x density four quark operators Qi by bosonized chiral 
fields from the chiral Lagarangian, then calculate loop contributions by using the functional 
cut-off momentum scheme. Finally, one matchs the two results obtained from QCD and 
ChPT with functional cut-off momentum by requiring scale independence of the physical 
results. 

The general form of the chiral Lagrangian can be expressed in terms of the expansions 
of the momentum p and quark mass to the energy scale A x . Here we only use the chiral 
Lagrangian which is relevant to the K — > tttt decays (for the most general one, see ref. JO] ) 
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C eff = J —{ tr(D^D"U) + — ^(lnt/* - \nU) 2 + r tr(MU ] + UM ] ) 

+r||tr (D^D»U(M ] U + f/ f .M) (5) 

x 

+r 2 ||tr (.M^A^C/ + A^^jM^) + r 2 ^tr(M^M) } 

x x 

with DfJJ = d^U — ir^U + iUl^, and M. = diag(m u , m^, m s ). and are left- and right- 
handed gauge fields, respectively. The unitary matrix U is a non-linear representation of 
the pseudoscalar meson nonet given as U = e m / f with II = n a \ a and tr(X a X b ) = 25 ab . Here 
we keep the leading terms at large N c limit except the anomaly term which arises from the 
order of 1/N C . Note that in order to make clear for two independent expansions, namely 
1/N C expansion characterized by p 2 / A 2 in the large N c limit, and the momentum expansion 
described by p 2 /A 2 c , we have introduced a scaling factor A x ~ 1 GeV and redefined the low 
energy coupling constants Lj introduced in ref. JL^] via Lj = x% / 2 /4A^ and Hj = Kj / 2 /4A x , 
so that the coupling constants x% = 3, 5, 8) and A x are constants in the large N c limit 
and the whole Lagrangian is multiplied by f 2 and is of order N c except the U(l) anomalous 
term. This is because when applying the large N c counting rules to all terms in Lagrangian, 
it shows that Li = 0(N C ) (i ^ 7) and Hj = 0(N C ). Here L-j is at order of N 2 . Numerically, 
one sees that Xi — 0(1) for A x = 1 GeV. 

As the first step, we represent the quark currents and densities by the chiral fields 

^ s = -*2 {u ' rTU 



X5_ 

2A X 



(d»tfM - M ] d»U + d^WUM^U - U^MU^d^u) ^ , (6) 

q jR q iL = - = - r £ + ^d^&*UlP + 2r^tfMtf + r^M^ (7) 

Similarly one can obtain the right-handed currents and densities. With these definitions, 
all the current x current and density x density four quark operators can be reexpressed in 
terms of the chiral fields, we may call such chiral representations of four quark operators Qi 
as chiral operators denoted by Qf correspondingly. At p 2 order, Qf can be written as the 
following form 



Q\ + H.c. = -f tr {XelPdfJJ) tr [X^lPd^U) + 0(1/ A 2 /) , 
Ql + H.c. = -f tr (wftdpUXMrfdPU) + 0(1/A 2 ) , 
qx + Ha = _f tr (x^d^u) tr (U^d^u) + 0(1/A 2 ) , 
qx + Kc _ = _f tr (x 6 d ll rfd»U) + 0(1/ A 2 X ) , 

qx + Kc _ = _f tr farfdpU) tr (u&>rf) + 0(1/A 2 ) , (8) 



QX + H.c. = +f (^jtr (X 6 d^d^u) + 0(1/Aj) 



qx + Ha = _ l _ Q x _ 3 /4 tr (x^d^u) tr (X^Ud^) + 0(1/ A 
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Q x s + H.c. = -\Ql + fr 2 \tr (X 6 U^^U 



2^° J 4 
.3 Xs 



+ / 4 r 2 -||trA 6 (u ] \ {1) Ud^d»U + d^U ] d»UU ] \ {1) U 

x 

+ fr 2 ~jj2rtr\ 6 (u ] 'X W U 'M^U + U ] ' MU ] X w u) + 0(1/Aj). 

x 



with the matrix A^^diag. (1,0,0). Thus loop contributions of the chiral operators Q} can 
be systematically calculated by using ChPT with functional cut-off momentum. 

For K — > 7T7T decay amplitudes and direct CP- violating parameter e'/e, the most im- 
portant chiral operators are Qf, Q%, Qq and Qg. In fact, the chiral operators and Q$ 
decouples from the loop evaluations at the p 2 order |§, i.e. 

Qt = QI = o (9) 

which can explicitly be seen from the above chiral representations due to the traceless factor 
tr = when ignoring the singlet U(l) nonet term which is irrelevant to the Kaon 

decays. Here Ud^W = A*\ a may be regarded as a pure gauge. This feature may also 
be understood as the fact that Q3 and Q5 operators are generated from strong penguin 
diagrams and suppressed by 1/N C factor in comparison with the operators Q4 and Qq, thus 
in the large N c limit, they decouple automatically. As a consequence, it implies that at the 
lowest order of p 2 , we arrive at two additional algebraic chiral relations 

qx = qx_qx = _ /4 tr (x^d^u) + 0(1 /A 2 X ) . (10) 

and 

Ql = - (^) {Qi - Qi) = (jr) f 4 ^ (Xedptftru) (11) 

Notice that the mass parameter r is at the same order of the energy scale A x , and X5 i s a t 
order of unit, thus the leading non-zero contribution of Q§ is at the same order of Q\ and 

QI 

The above algebraic chiral relations were first derived in ref. [0, they have also been 
checked from an explicit calculation up to the chiral one-loop level by using the usual cut-off 
regularization |19| . If naively identifying the cut-off momentum M to the QCD running scale 
/j,, the above algebraic chiral Wu-relations, as commented by Buras, Jamin and Lautenbacher 



24], seem to hold only at one point when matching to QCD. Thus two questions have been 



arised: 

• From which energy scale and up to which order of chiral loop corrections the algebraic 
chiral Wu-relations hold; 

• How can the algebraic chiral Wu-relations survive when matching the ultraviolet cut- 
off momentum of ChPT to the infrared cut-off momentum of perturbatve QCD. 

Let us briefly issue the first question and leave the second question to next section. The 
answer to the first question is manifest, the algebraic chiral Wu-relations hold starting from 



7 



the energy scale where the low energy dynamics of QCD is considered to be described by 
the ChPT with functional cut-off momentum. They even survive when we include chiral- 
loop corrections generated from the lowest p 2 order terms as Wu-relations are the algebraic 
chiral operator relations, they should not be modified by the chiral loops of the strong 
interactions. Note that the coupling constants must also be replaced by the renormalized 
ones at the same order oil/N c or 1/A^. The reason is simple as the non-trival structures of 
Qg and (Q* — Q\) at the order of p 2 and p 2 /N c (or p 2 / K 2 F ) are unique. The order of p 4 terms 
are suppressed by the factors m 2 K /A 2 c and m 2 /A 2 . Therefore, in the chiral limit, namely 
m 2 K , m 2 <C A 2 the above algebraic chiral operator relations should hold up to the order of 
p 2 /N c and p A . This may be understood in an analogous way to QCD, where relations of the 
quark operators survive from all order of QCD corrections. The reason is simply due to that 
QCD is a renormalizable theory. For the ChPT, though it is an effective theory and not a 
renormalizable one in the usual sense, but it can be constructed to be a consistent theory 
order by order in the expansion of momentum and quark mass as well as 1/N C . Thus, to a 
given order of expansion, ChPT may be regarded as a renormalizable one in the more general 
sense ||25|| . Therefore, the algebraic chiral operator relations must survive, at least, up to 
the one-loop corrections, which has actually been checked from our explicit calculations. It 
was based on this observation, we came to our early conclusion that the direct CP- violating 
parameter e' / e can be large enough to be measured and its favorable numerical value is likely 
to be around e'/e ~ (10 — 30) x 10~ 4 PJlQf. As the matching to QCD was not completely 
considered, our previous results |9|Jl(| strongly depend on the strange quark mass. 



IV. MATCHING BETWEEN QCD AND CHPT 

Let us begin with the short-distance operator evolution from perturbative QCD. When 
the energy scale \i is high, mw > £t > rrib, there are eleven independent operators Qi 
(i = 1, • • • , 11). When the energy scale \i runs down to below the bottom quark mass m& 
and above the charm quark mass m c , i.e., m b > ji > m c , the operator Q u decouples and 
operator Q w is given by the linear combination Q w = —2Qi + 2Q 2 + Q3 — Qi- Once the 
energy scale /1 goes down to below m c but above the confining scale or the energy scale A x , 
i.e., m c > /i > A x , two operators Qg and Q4 become no longer independent and are given 
by the linear combination Qg = Q 2 + Qi and Q4 = Q3 + Q2 — Qi- Thus there are only 
seven independent operators below m c and above A x . In order to match to the long-distance 
evolution of the operators, let us present one-loop QCD corrections of the quark operators 
at the energy scale just above the energy scale A x 

2 

QiQki) = QM ~ 3^ ln(^f ) Q a (Ai) + 0(1/N C ) , (12) 



2 

Q 2 (^)=Q 2 (/i)-3^1n(^f)Q 1 ( yU ) 

l £ ln( f } - \ S H f } QM + 0{1/Nc) ' (13) 



S 



,2 „ ,,2 



" £ ln( 7^ g4(/i) " S ln( ^ } QM + 0{1/Nc) ' (14) 



,2 



Q 6 (/iQ) = g 6 (/i)-^in(^f)g 4 (/i) 

47T /i z 



,2 



+ [3(N C - 1/N C ) - 1]^ ln(^) g 6 (/i) + 0(l/iV c ) , (15) 



,2 



Qsfo) = QM + [3(JV C - l/JV e ) - 1]^ ln(^) QM , (16) 



and 



11 a s , ,^L n , v 2a s Jig, 



Qs(^) = QsM - T -p In(^) g 4 (/x) - - / ln(^) QM + 0(1/N C ) , (17) 

2 

g 5 (/i Q ) = QM + 3^ ln(^f ) g 6 (/i) + 0(1/N C ) , (18) 

47T /i z 

2 

g 7 (/i Q ) = g 7 (/i) + 3^ in(^fi) g 8 (/i) + o(i/n c ) . (19) 

From the above results, we come to the following observations: 

i) In the large N c limit, Q±, Q2, Q4, and Qq form a complete set of operators under QCD 
corrections. 

ii) The evolution of Q$ is independent of other operators and only caused by loop cor- 
rections of the density. 

iii) The operator Q 3 is given by the linear combination Q 3 = g 4 — (g 2 — Qi). The 
operator Q 5 is driven by the operator Q 6 , and the operator Q 7 is driven by the operator Q s . 

When the energy scale \i approaches to the confining scale, or /1 < A x ~ ~ lGeV, 
as we have discussed in the above sections, long-distance effects have to be considered. The 
evolution of the operators Qi{^) when /1 < A x is supposed to be carried out by the one of 
the chiral operators Q*(M(fi)) in the framework of the functional cut-off ChPT truncated 
to the pseudoscalars. To be treated at the same approximations made in the short-distance 
operator evolution of QCD, we should only keep the leading terms (i.e., quadratic terms 
of functional cut-off momentum) and take the chiral limit, i.e., m 2 K ,m1 << A 2 F . In such a 
leading 1/N C approximation and chiral limit, we find that the evolution of the operators Q\ 
and Q\ can be simply given by the following forms when the functional cut-off momentum 
runs from M(/i) down to M(fjf) 

QM - g?(MGu)) = Qj(MOiO) - 2(M2(/i) A ; MV)) Qm^')) , (20) 

QM - QiiMiti) = QmM) _ 2 ( M2 M - M2 (^)) gx (W)) 



A 



mm - MM) 



+ ,2 ^ (Q2 ~ Q X i)(M(»')) , (21) 



A-F 



where A^ = AnF = 1.16 GeV with F the renormalized one of /. Notice that the operators 
Qz and Q^ decouple from the evolution, namely Q% = and Q% = 0. The results for the 
operators Q* (i = 4, 6, 8) can be written as follows 
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QM - Qi{M{^)) = (Q x 2 - Q?)(m(aO) , (22) 

QM - Q*(//, M(//)) = + 3(iV c - 1/iVc)^ ln(^)j Q£(// x , M(//)) , (23) 

Q 8 (//) - qftdi, M{y)) + 3(N C - l/JV c )g ln(j^) j Q* (/i x , M(//)) . (24) 

where the explicit //-dependence of the operators Qq(/i, M(fi)) and Qg(n, M ((/,)) arise from 
the running quark mass and behavors like l/(m s (/i) + m(/t)) 2 . Notice that the independent 
operators are reduced once more in the long-distance operator evolution when \x < A x due 
to the algebraic chiral operator relations. Let us now compare and match the loop results 
evaluated from QCD with the ones from the ChPT with functional cut-off momentum at the 
energy scale A x . Substituting eqs.(20) and (21) into eq.(12), keeping the leading 1/N C terms, 
we obtain, from the requirement of //-independence in the large N c limit, i.e., -§-Qi(^q) = 0, 
the first matching condition, 



d_ f 2M 2 (jj) \ _ 3a, 



^R^ =S, (25) 



which can in general be obtained by requiring the matching between the anomalous dimen- 
sions of quark operators Qi(fi) in QCD and the ones of the corresponding chiral operators 
Qf(M(fi)) in ChPT in the large N c limit, i.e., 

if eson = ^QUM(ri) = i? uark = ^Qi(^) ■ (26) 

Analogously, substituting eqs.(20)-(23) into eq.(13), keeping the leading l/N c terms and 
adopting the above first matching condition, we arrive at the second matching condition 

C»(M x ,M( A i)) = ~(^f_Qf)(M( A i)), //<A X (27) 

Note that such a matching condition holds for whole energy scale \i < A x . At the special 
point M(/i) = 0, it covers the condition first presented in ref. ^(J. In fact, the above two 
matching conditions may simply seen by comparing eqs.(20) and (21) with eqs.(12) and 
(13). On the other hand, from the chiral representation of operators and their chiral loop 
corrections, we have the following chiral relation in the leading 1/N C approximation and 
chiral limit 

Q*Gu x , M(ji)) ~ ("^) (Q2 ~ Q?)MaO) , A* < A x (28) 

R x = R(ji ~ fi x ) ~ m 2 /m(/i x ) ~ 2m 2 K /(m s + m)(ji x ) , (29) 

where we have simply replaced the coupling constants %5 and r by the corresponding renor- 
malized ones x\ an d -R(a0 as their loop corrections are at the subleading order. When 
combining the second matching condition, it allows us to fix the strange quark mass 



Kxl _ ii 

A x - 2 



m s (// x ) ~ 196MeV , (30) 
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Here we have used the result A x = 1.03^X5 GeV which is fixed from the ratio of the kaon 
and pion decay constants. 

The first matching condition can be rewritten as follows after integration 

^ = !^ + ^ ln( 4), (31) 
A% A 2 F 4tt Vo 

where /io and M = M(/i = /i ) are two integral constants. It is seen that the /i-dependence 
of the functional cut-off momentum M(/i) is now logarithmic. Noticing the approximation 
ln(/i 2 //iQ) ~ H 2 /Ho — 1 when /i 2 ~ /1q, namely the functional cut-off momentum M(/i) 
is approximately proportional to ji when /i runs down and approaches to the low energy 
scale /i which is expected to be slightly above the QCD scale Aq C d, we then have Mq = 
(3a s (fi )/8ir)A 2 F . Thus the /i-dependence of the functional cut-off momentum M(/i) can be 
written as 

Al ~ 4vT + 4vT in V§ j ' {62) 

which shows that after imposing the matching condition for the anomalous dimensions be- 
tween quark operators Qi(fi) in QCD and the corresponding chiral operators Q*(M(/i)) in 
the ChPT with functional cut-off momentum, the dimensionless ratio M 2 / A 2 F is only related 
to the strong coupling constant a s and becomes scheme- independent, which implies that the 
long-distance operator evolution in ChPT with functional cut-off momentum can be carried 
out by using any approach. For instance, with and without separating factorized and non- 
factorized contributions, should obtain the same results after appropriately considering the 
matching between QCD and ChPT with functional cut-off momentum. 

In general, we have /io > Aqcd- To fix the value of /io, we use M ~ /io- Thus jiQ (or 
a s (/io) ) is determined via 

/i ^ A F ^3a s (fi )/8n (33) 

Using the definition a s (/i) = 67r/[(33 — 2nf) ln(/i 2 / Aq CD )] with rif = 3, the initial low energy 
scale /i is found, for Aq C d = 325 ± 80 MeV, to be 

/i ^ 435 ± 70MeV or a s (^)/2n ~ 0.19l{J;{jjj . (34) 

With such an initial value of /i , the functional cut-off momentum M(/i) at fj, — A x yields 
the following corresponding value 

M x = M(n = A X ~ lGeV) ~ 0.71±%£GeV . (35) 

which provides the possible allowed range of the energy scale where the ChPT with functional 
cut-off momentum can be used to describe the low energy behavior of QCD at large iV c limit. 



V. EVOLUTION OF LONG-DISTANCE CHIRAL OPERATORS 

From the above analyzes, the AS* = 1 low energy (/i < A x ) effective Hamiltonian for 
calculating K — > nn decay amplitudes may be written as 
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V% , | . 2 . J (>> 



ntff^^K £ Qm(ri) , (^<a x ) (36) 



we may now adopt the matching conditions and algebraic chiral operator relations to in- 
vestigate the evolution of the chiral operators Qf(M(A x )). The first matching condition 
enables us to sum over all the leading terms via renormalization group equation down to 
the energy scale /x , and the second matching condition together with the algebraic chiral 
operator relations allows us to evaluate the penguin operators Q*(M) and Qq(M) from the 
operators Qf{M) and Q%(M). So that the operators Q\(M) and Q%(M) form a complete 
set for the operator evolution below the energy scale fi ~ A x ~ 1 GeV, or correspondingly 
below the functional cut-off momentum M(fi ~ A x ) ~ 0.7ll"o;i2 GeV for Aqcd = 325 ± 80 
MeV. It is convenient to choose a new operator basis Q±(M(/i)) = Qfe^M^ij)) ± Q^{M{ix)). 
The anomalous dimension matrix for the basis (Q*, Q\) is found to be 

a s f -9/2 \ f . 

7 = Y (-3/2 3 J ^ 

Following the standard procedure of the renormalization group evolution with the initial 
conditions for the Wilson coefficient functions: c_(A x ) = C2(A X ) — ci(A x ) and c + (A x ) = 
c 2 (A x ) + Ci(A x ) , we find in the leading logarithmic approximation that 

g^(M(A x ))=r ?x 1 /2gx (/io) > (38) 
Q^(M(A X )) = v ^Ql fa) + I fc V2 " vT) Q-(»o) , (39) 
with r\ x = a s (A x )/a s (/i ), and 

Q*(/i )=Q*(0) + ^^Q*(0), (40) 

Q$M = QW) - 3 -^± q$(o) + (o) , (4i) 

In the above analyzes, we have taken the chiral limit m 2 K ,m^ <C K 2 F . From the chiral 
one loop results, the finite meson mass contributions can be approximately included by 
modifying the above results into the following form 

Q*(M(A X )) =r,-^r,_(M x )Q x -M , (42) 

Q\{M{K X )) = ^ m (M x ) Qiino) + I (^ x 1/2 - vT) V2(M x )Q*to) , (43) 

with 

*W^l + =^b.(l + ^). (44) 
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Numerically, we use m ~ 300MeV, m K = 0.495 GeV and m* = 0.137GeV. When the QCD 
scale takes the value Aqcd — 325 ± 80 MeV with the corresponding low energy cut-off 
momentum /i — 435 ± 70 MeV, we have 

QHM(AJ) = (3.171^) Q*(0) = g 4 x (M(A x )) , (45) 
Ql(M(A x )) = (0.55^;") Ql(0) + (0.81^) C£(0) , (46) 

Q^, M(A X )) = -Hgx (M(Ax)) = -(17.4411^) Q* (0) , (47) 

33 A 2 

Qii^ M(A X )) = £ * (Q\ + QPL)(0) = 19.18 {Q\ + C£)(0) . (48) 

Which show that the isospin 1 = 2 amplitude A 2 is suppressed by a factor of about 2 as it only 
receives contributions from the operator g*(M(A x )), while the isospin 1 = amplitude A 
is enhanced by a large factor as it mainly gets contributions from the operator g* (M(A X )). 
On the other hand, the direct CP- violating parameter e'/e is expected to be large since the 
significant enhancement of Q\ (/i x , M(A X )) relative to Q\ (/z x , M(A X )) is seen to be resulted 
from the algebraic chiral operator relation and matching condition. We are going to present 
our numerical predictions for the isospin amplitudes and the direct CP-violating parameter 
e'/e in the next section. 

VI. PREDICTIONS FOR e'/e AND AI = 1/2 RULE 

We are now in the position to calculate the K — > nn decay amplitudes Aj with isospin I 
A I cos6 I = (nn\H?fj 1 \K) = ^\ u £ c i (A x ) J Re(g^(M(A x ))) / (49) 

i=l,2,4,6,8 

The CP-conserving amplitudes are given by 

ReA cos5 = ^Re\ u Zl (A x )Re(Qf(M(A x ))) 

V 2 i=l,2,4,6,8 

~ ^ Re\ u [±z-(A x )Re{Q*(M(A x )))o + \z + {A x )Re{Q\{M{A x ))), 
+ z,(A x )Re(Q*(M(A x ))) + z Q {A x )Re{Ql{M{A x ))) ] , (50) 
ReA 2 cos5 2 = ^Re\ u ]T Zl {A x )Re{Q}{M{A x ))) 2 (51) 

V2 i=l,2,8 

~ ^ Re\ u [\z_{A x )Re(Q*{M{A x ))) 2 + \z + {A x )Re{Q\{M{A x ))) 2 \ , 
and the CP-violating amplitudes are dominated by (Qq(M(A x )) and (Qg(M(A x )) 2 
/mA cos(5o = /mA t ]T Vi(K)MQ? (M(A X ))) 

v2 i=l,2,4,6,8 

~ 7mA t [t/ 6 (A x ) J Re(g^(M(A x ))) ] , (52) 
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ImA 2 cos5 2 = -^ImX t £ yi (A x )Re(Qf (M(A X ))) 2 

v2 i=l,2,8 

~ ^ JmA t [y 8 (A x )i?e(Q^(M(A x ))) 2 ] . (53) 



From the definition of direct CP- violating parameter e'/e 

e' uj ( ImAz ImAo 



e y/2\e\ \ReA 2 ReA . 
with uj = ReA 2 / ReAo = 1/22.2, we arrive at the following general expression 

e' G F uj 



(54) 



2 \e\ReA 



ImX t (h - h 2 /uj) (55) 



Here ho and h 2 are given by the isospin 1 = and 1 = 2 hadronic matrix elements of relevant 
operators 

ho = (cos^o)- 1 £ yi (A x )Re(Q? (M(A X ))) (1 - Q IB ) 

i=l,2,4,6,8 

~ (cos5 )- 1 y 6 (A x )^e(Q^ (M(A X ))) (1 - Q^r) , (56) 
h 2 = (cos6 2 )- 1 ]T y l (A x )Re(Qf(M(A x ))) 2 

i=l,2,8 

~ (cos5 2 )- 1 2 / 8 (A x ) J Re(g^ (M(A X ))) 2 , (57) 

where we have taken into account the possible isospin breaking effect Qib, its previously 
estimated value was Qib — 0.25 ± 0.1 f26|| . The most recent refined calculation in |27| gives 



a smaller value Qib — 0.16±0.03 but with a large error [|28||. In our present numerical calcu- 



lations, we use Qib — 0.16. The CKM factors Re\ u and Im\ t are given in the Wolfenstein 
parameterization as follows 

i?eA u = Re(V: s V ud ) = X , ImX t = Im(V*V ta ) = A 2 X 5 V (58) 

To evaluate the numerical results, we are going to take the following reliable values 
for all relevant parameters. For the involved energy scales, we have: Aqcd — 325 ± 80 
MeV, no = 435 ± 70MeV, A x = l.OGeV and A F = 1.16 GeV. For the Wilson coefficient 
functions, we only use the leading order results at one-loop level for a consistent analysis 
since the chiral operators have only been evaluated up to the leading order at the chiral 
one-loop level, namely at the order of 1/N C ~ M 2 /A F ~ a s in the large iV c approach. Their 
values can be read following refs. |13],[T4]]. The numerical values at fi = A x are regarded 
as the 'initial conditions' for the chiral operator evolution and read for Aqcd — 325 ± 80 
MeV : z_(A x ) = (z 2 - Zl )(A x ) = 2.18118$?, z + (A x ) = (z 2 + Zl )(A x ) = 0.685 =F 0.029, 
Z4 (A X ) = -(0.012±0.003)and^ 6 (A x ) = -(0.013±0.003), as well asy 6 (A x ) = - (o.H3l^) 
and ys(A x )/a = 0.158t[J;o33. The hadronic matrix elements of chiral operators at cut- 
off momentum M = take their values at the tree-level: (Q*(0)) = 36.9 x 10 6 MeV 3 , 
(Ql(0)) = 12.3 x 10 6 MeV 3 , (Q$(0)) 2 = 34.8 x 10 6 MeV 3 and (Q£(A x ,0)) 2 = 328.8 x 10 6 
MeV 3 . For the CKM matrix elements, there remain big uncertainties arising from the single 
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CP-violating phase, two matrix elements V u \, and V c b, or the corresponding Wolfenstein 
parameters 77, p and A. For a numerical estimate, we take Re\ u = 0.22 and Im\ t = 1.2 x 10~ 4 
as the central values. With these input values, we obtain, in the leading 1/N C approximation, 
the isospin amplitudes 

ReA = (2.56l°;$) x 10" 4 (cosc^ 1 MeV = (3.10^ 4 ) x 10" 4 MeV , (59) 
ReA 2 = (0.12 t 0.02) x 10~ 4 (cos5 2 ) _1 MeV = (0.12 t 0.02) x 10~ 4 MeV , (60) 

which agree well with the experimental data: ReAo = 3.33 x 10~ 4 MeV and ReA 2 = 
0.15 x 10 -4 MeV . Here we have used the final state interaction phases, 5q = (34.2 ± 2.2)° 



and 5 2 = (—6.9 ± 0.2)° |j30| . Consistently, the direct CP- violating parameter e'/e is found, 



in the leading 1/N C approximation with Aqcd = 325 ± 80 MeV, to be 

e L = (23.6l^ 8 4 ) x 10~ 4 (61) 

for the central value of Im\ t = 1.2 x 10~ 4 resulted from fitting the indirect CP-violating 
parameter e, \V u b\, B° — B° and B° s — B° s mixings. It is of interest to note that this central 
value also agrees with the one predicted from ten useful relations among the quark masses 
and mixing angles obtained in the SUSY GUT model ||31|| . When considering the possible 
allowed range for the CKM matrix elements extracted from fitting the present experimental 
data, we have 



(13-8±I;|) x 10- 4 < t = (23.6^ 2 8 4 ) x 10" 4 ( J™^ ) < (33.6^^) x 10" 4 (62) 



for the possible allowed range 0.7 x 10~ 4 < Im\ t < 1.7 x 10 -4 . It is noticed that the present 
new predictions for the isospin amplitudes and direct CP- violating parameter e'/e further 
confirm our early conclusions PJlOfl. Our new predictions are consistent with the world 
average || 

Re{e'/e) = (19.2 ± 2.4) x 10~ 4 (World Average 2000) (63) 

which is obtained by taking into account the results from four collaboration groups. They 
contain two published results reported early by NA31 collaboration and E731 collaboration 

Re(e'/e) = (23 ± 7) x 10" 4 , (1993 NAZI) [4] ; 

Re(e'/e) = (7.4 ± 5.9) x 10" 4 (1993 £731) § (64) 

and the recent new results reported by the KTeV collaboration at Fermilab and the NA48 
collaboration at CERN: 

Re(e'/e) = (28.0 ± 3.0 ± 2.8) x 10~ 4 (1999 KTeV) [6] ; 
Re(e'/e) = (18.5 ± 4.5 ± 5.8) x 10~ 4 , (1999 NAA8) § ; 
Reie'je) = (14.4 ±4.3) x 10" 4 (2000A^A48) § , (65) 

Before drawing our conclusions, we would like to address the following points: 
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1. The main uncertainties for the predictions arise from the QCD scale Aqcd (or the low 
energy scale /xq) and the combined CKM factor Im\ t . Nevertheless, the uncertainties from 
the energy scale may be reduced from comparing the predicted isospin amplitudes A and 
A 2 with the well measured ones. As a consequence, it is seen from eqs.(59) and (60) that 
the results corresponding to the large values of Aqcd > 325 MeV appear not so favorable. 

2. From the above point of view, it is seen that from the isospin amplitude A 2 , the ratio 
e'/e favors the low values 

pi 

- = (16 ± 7) x 1(T 4 (66) 
while from the isospin amplitude Aq, it favors the high values 

- = (24 ±10) x 10~ 4 (67) 

From the ratio of the two amplitudes ReA / ReA 2 , i.e., the A J = 1/2 rule, the ratio e'/e 
favors the middle values 

- = (20 ± 9) x 10~ 4 (68) 

which is consistent with the most recent results reported by the NA48 collaboration at CERN 
IK?] and the KTeV collaboration at Fermilab |J. In fact, it is very close to the average value 
from NA48 and KTeV. While the central values from two experimental groups differ from 
each other at 3.5-cr level. 

3. The above results are renormalization scheme independent as the consistent matching 
between QCD and ChPT considered in the present paper is at the leading one-loop order 
of 1/N C ~ a s ~ 1/Afi around the scale A x . The renormalization scheme dependence arises 
from the next-to-leading order of perturbative QCD, which could become substantial for 
some of the Wilson coefficient functions when the renormalization scale /x runs down to 
around the scale A x = lGeV. In our present approach, the scheme for the long-distance 
evolution is fixed by the ChPT with functional cut-off momentum. For matching to this 
scheme, it is useful to introduce a scheme independent basis for the perturbative QCD 
calculation of short- distance physics. Then applying our above procedure to find out the 
matching conditions at the next-to- leading order l/Af? ~ a 2 s ~ l/A 4 ?. To work out the 
scheme independent basis, it may be helpful to adopt the method discussed in ref. [0 and 



use the cut-off momentum basis. The study of scheme independent basis in perturbative 
QCD is beyond the purposes of the present paper. Some effort is being made |33| though it 
is not yet full understood. 



VII. CONCLUSIONS 

We have simultaneously analyzed the direct CP- violating parameter e'/e and AI =1/2 
rule in kaon decays by considering a consistent matching scheme between QCD and ChPT. 
Our main points may be summarized as follows: i) Starting from the chiral Lagrangian 
obtained in terms of the momentum and quark mass expansion with low energy coupling 
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constants given by the leading terms of the 1/N C expansion. The N c behavior in the con- 
cerned chiral Lagrangian has explicitly been characterized by the scale factor ~ N c ; ii) 
The four quark operators for weak kaon decays at low energies have been assumed to be 
represented by the chiral operators in the large N c limit. It has also been shown that there 
is simplification in ChPT which leads to useful algebraic chiral operator relations. Those 
relations survive even when including loop corrections; iii) We have adopted the usual cut- 
off regularization scheme Jill for all the diagrams with a single cut-off momentum for a 



systematical analysis, and given up the scheme of separating factorized and non-factorized 
parts with two cut-off scales; iv) The cut-off momentum M has been considered to be the 
function of the QCD running scale n, i.e., M = M(/x), instead of naively identifying it to 
the perturbative QCD running scale /i. The form of the functional cut-off momentum M(/x) 
has been determined via the matching requirement, so that the chiral loop results become 
scheme independent. As a consequence, two useful matching conditions have been obtained, 
which has allowed us to evaluate the long-distance chiral operators and sum over the lead- 
ing non-perturbative contributions. Of particular, the AI = 1/2 rule can consistently be 
understood and the resulting direct CP- violating parameter e' / e become large enough to be 
measured, which also confirms our early conclusions PJIOfl. Taking into account the simul- 
taneous consistent analysis for the isospin amplitudes A and A 2 , the numerical result for 
the ratio e'/e is in favor of the values 

£ > 

- = (20 ± 9) x 1(T 4 (69) 

£ 

which may be regarded as the favorable prediction in our present analyzes. The prediction is 
also consistent, within the theoretical and experimental uncertainties, with the present data 
HUHZtl- Finally, we would like to remark that we have neglected in our present analyzes the 
subleading contributions, their effects are in general small and will be investigated elsewhere 
in detail. 
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